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Abstract 



In this paper we propose a generalization of the Gilkey-de Witt heat kernel expan- 
sion, designed to provide us with a precise estimation of the heat trace of non-negative 
Schrodinger type differential operators with non-trivial kernel over all the domain of its 
"inverse temperature" variable j3. We apply this modified approach to compute effectively 
£f) ' the one-loop kink mass shift for some models whose kink fluctuation operator spectrum 

is unknown and the only alternative to estimate this magnitude is the use of the heat 
kernel expansion techniques. 

PACS: 11.15.Kc; 11.27.+d; ll.10.Gh 

<D 

£j- 1 Introduction 



In the mid sixties outstanding developments concerning the high-temperature expansion of 



the kernel of the generalized heat equation associated with a differential operator of the 
Laplace or Dirac type took place in different fields from Mathematics [H[2] and Physics [3]. In 
particular, Gilkey unveiled the meromorphic structure of the spectral zeta function, showing 
■ that the residua at the poles of this function are determined from the Seeley coefficients of the 

heat kernel expansion and describe certain topological invariants/characteristic classes. On 
the other hand De Witt used the heat kernel expansion to deal with quantum fields on curved 
backgrounds. Thereafter the heat kernel/zeta function methods have become an important 
tool in Quantum Field Theory [H [6j [7] . 

The computation of the one-loop kink mass correction is a topic enclosed in this frame- 
work. Technically this magnitude quantifies the contribution to the mass of the quantum 
fluctuations over the classical kink solution measured with respect to the background of vac- 
uum quantum fluctuations, the kink Casimir energy. A mass renormalization must also be 
implemented in the previous scheme. All of this can be expressed in terms of the spectral 
zeta functions of the kink and vacuum fluctuation operators, which in turn can be written as 
an integral of its heat traces over the interval (3 G [0, oo) by applying the Mellin transform. 
During the last decade we (and our colleagues) succeeded in applying the Gilkey-de Witt 
(GDW) heat kernel expansion in (1 + l)-dimensional scalar field theoretical models to ex- 
press the kink mass shift as a truncated series in the Seeley coefficients of the kink fluctuation 
operator heat trace, see [3 [9]. In some applications -heat kernel proofs of the index theorem, 
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computation of anomalies in QFT- the relevant information is encoded in the behavior of 
the heat trace for small values of the "inverse temperature" variable /?. The kink mass shift, 
however, demands a precise estimation of the heat trace over all the domain of the variable, 
P € [0,oo). We recall that the small (3 regime captures the higher eigenvalues therefore 
determining the ultraviolet behavior while the large f3 range is contrarily dominated by the 
lower eigenvalues and it is important in the infrared domain [TO]. The previous assertion is 
of no consequence if we dealt with strictly positive operators because in this case the GDW 
approach is well established; however it is crucial for the kink fluctuation operator, which 
always comprises a zero mode and consequently the GDW expansion does not reproduce the 
asymptotic behavior of the heat trace adequately. The usual manoeuvre to solve partially 
this problem is to truncate the Mellin transform in the range where there exists a good fitting 
between the GDW expansion and the heat trace. This obviously leads to a loss of precision 
in the final computations. 

In this paper we introduce a different alternative; we modified the GDW route to accom- 
modate the role of the zero modes in the theory. This is described in Section §.2. In Section 
§.3 we shall apply this modified GDW heat kernel expansion to obtain a precise computation 
of the one-loop kink mass shift for a (1+1) dimensional scalar field theoretical model, where 
other more direct methods (DHN formula, see [H]), are unapproachable because the kink 
fluctuation operator spectral information is unknown. 



2 The modified heat kernel asymptotic expansion 
2.1 Heat kernel asymptotic expansion and zero modes 

In this section we shall address the generalization of the GDW heat kernel asymptotic ex- 
pansion for non-negative second order differential operator of the form 

d 2 

K = -—r + v 2 + V(x) xeR , (1) 

dx A 

where V(x) is a real function whose asymptotic behavior complies with 

lim V{x) = . 

x— >-±oo 

In this case the iT-spectrum 

Spec(iT) = {ul = 0} U {uj 2 n } n =i,...,i U {k 2 + v 2 } keR 

embraces a zero mode fo(x), £ bound states f n (x) with non-negative eigenvalues and scat- 
tering states fk(x) emerging on the threshold value v 2 . The computation of the heat integral 
kernel by means of its definition 

1 

K K (x,y;P) = f^y)fo(x) + ^2f*(y)f n (x)e-^ + / dk f* k (y) f k (x) (2) 

n=l ^ 

is not possible for the most of the cases because we generally lack the spectral information 
of the JT-operator. An alternative way to estimate Kx{x,y; /3) is to exploit the fact that ([2]) 
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is a solution of the parabolic (heat) equation 



dj3 dx 2 
with the contour conditions 



9 d2 2 \ t ( ^ 



K K (x,y;P) = , (3) 



lim K K (x,y;P) = fo(y)fo{x) , lim K K (x, y; fi) = S(x - y) , (4) 

dictated by the definition ([2|). Therefore the zero mode rules the asymptotic behavior of 
Kk(x, y; j3). This constitutes a variation from the original GDW approach, where the second 
hand side of the first relation in @ vanishes because of the positiveness assumption upon 
the -fT-spectrum assumed in that framework. 

The potential wells of the -RT-operator (JTj) and the Helmholtz operator 

d 2 

Ko = ~dx 2+y2 , xeR , (5) 

have the same asymptotic behavior. The spectrum of this latter operator Spec(-Ko) = {k 2 + 
v 2 }kt=M. only comprises scattering states = -^=e tkx . From this the exact i^o-heat kernel 

is 

K KQ {x,y^) = -^=e-^ 2 e- {J ^ L , (6) 
whose contour conditions verify 

lim K Ko (x,y;f3) = , lim K Ko (x, y; /3) = 5(x - y) . (7) 

,3-»+oo p^fO 

The original GDW method analyzes the /f-heat kernel Kk(x, y, f3) by means of its deviation 
from the exact form ([6]) using the standard factorization 

K K (x,y;P)=K Ko (x,y;P)A{x,y;l3) with A{x,y;0) = l . (8) 

The relation ([8]) is, however, inconsistent in the present context because the two sides of ([8]) 
have different j3 — > oo asymptotic behaviors, see (HD and ([7]). The reason of this underlies the 
presence of the zero mode in the if-spectrum. Notice, however, that ([8]) can be used to study 
the heat kernel Kx(x,y; f3) in the small (3 range. This make possible, for instance, the heat 
kernel proofs of the index theorem, the computation of anomalies in QFT, etc. Nevertheless, 
other applications, such as the computation of the one-loop kink mass correction, demand the 
integration over the interval f3 G [0, oo) of the heat kernel trace. Consequently the standard 
factorization ([8]) must be modified to reproduce the behavior of this function in the large f3 
range. 



2.2 The modified Gilkey-de Witt heat kernel expansion 

For non-negative differential operators of the type ([1]) the factorization ([8|) must be replaced 
by the form 

K K (x,y-,(3)=K Ko (x,y-,P)C(x,y-,P)+g(P)e-^fZ(y)f (x) , (9) 
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in order to adapt the GDW formalism to the existence of zero modes fo(x). The expression 
@ complies with the contour conditions provided that 



lim C(x, y: 0) = 1 



lim g(fi) = 1 , lim g(J3) = . 
0— >oo /?— >0 



Plugging the ansatz ([9]) into (J3]) leads to the "transfer" equation for C{x,y;j3): 

d 2 







5 x — y d 
d(3 (3 dx dx 2 



+ V(x) C(x,y;p) + 



+ v / 4^e^ / *(y) 



(10) 

(11) 

(12) 



The terms specified in the line (|12p concern to the role of the zero mode in this new scheme 
and complement the standard GDW "transfer" equation (written in the line (jlip). The PDE 
(|11|) is traditionally solved by means of a power series expansion: 



C{x, y; /3) = ^2 °n(x, y) P n with c (x, y) = 1 



(13) 



n=0 



where the constraint co(x,y) = 1 is derived from the first condition in (|10p . Plugging the 
expression (fT3j) into (jXTJ) the relation 

d 2 c n (x,y) 



n=0 



(n+ l)c n+ i(x,y) 



dx 2 



holds. A significant simplification of ()14[) is obtained by setting 

=> g(P)=Er£(vy/p) , 



dp 



(14) 



(15) 



which complies with the asymptotic conditions (|10p imposed on the function g(f3). Now 
the GDW procedure can be followed by implementing the asymptotic expansion of the error 
function 



Er fz = ^=e- z2 y 



n=0 



(2n + l)! 



■z 



2n+l 



in the expression (|14|) . This leads to the recurrence relations for the coefficients c n (x,y): 

n t . i\ t \ d 2 c n {x,y) dc n+l {x,y) . . . 
= (n + l)c n+ i(x,y) ^ H^-y) HI h F(x)c„(x, y) + (16) 



+2vfS(y)f(x)S 0n + f5(y)f(x] 



2^+1^271+1 



<9x 

d/o(x)2"+V-+i 
(2» + l)!! +( g -^o(y) dx (2n + 1)!! ( 17 ) 

where again the terms displayed in (|17p describe the new contributions to the standard 
expression (|16|) . From the factorization Q, the power series expansion (|13p and the choice 
()15p of g((3), the heat kernel can be written as the series expansion: 



K K (x,y;f3) = K Ko (x,y;f3) £ c n (x, y) f3 n + Erf ^e'^ f^y)f 



(18) 



n=0 



where we recall that fo(x) is the zero mode of the .ff -operator and the coefficients c n (x,y) 
are computed by using the recurrence relation (|16p starting from co(x,y) = 1. 
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2.3 The modified GDW heat trace expansion 

From the diagonal of the heat kernel 

-fSv 2 00 

K K (x,x,0) = \im K K (x,y,l3) = -== V c n (x, x) f3 n + Erf(^)/ *(x)/ ( 

y->x \/4vrp ' 

v n=0 

we can derive the spectral if-heat trace = Tr^2 e - ^ as 



(19) 



h K W) 



dx Kk{x, x; j3) . 



The Seeley densities c n (x,x) introduced in (fP9"j) are defined as 

c n (x,x) = lim c n (x,y) . 

y->x 



(20) 



(21) 



From (|16|) a recursive relation can be constructed for these new coefficients. To do this it is 
necessary to deal with the following subtlety: the operations of taking the y — > x limit and 
the derivatives with respect to x in the formula (|16|) do not commute. To cope with this 
problem we introduce the new auxiliary coefficients: 



whose first coefficients ^Cq(x), k E N are 



^C (x) = lim dkc ° - * k0 



dx k 



(22) 



(23) 



Taking the fc-th derivative of (|16|) with respect to x and later passing to the y — > x limit, the 
rest of the coefficients ^C n (x) verifies: 



{k) C n (x) 



1 



n + k 



IH,c -'W-E( J ')| M u.) 



-2^)^^-/0(^^(1+2*) 



which must be used in a recursive way until the Seeley densities in ([19 

c n (x,x) = ^C n (x) 



(24) 



(25) 



are determined. The first three Seeley densities c n (x, x) derived from f)24|) for the i^-operator 
are listed: 



co(x,x) = ^C {x) 
Cl (x,x) = C°)Ci(s) 

c 2 {x,x) = ^C 2 {x) 



1 



-V(x) 



±vf Q {x) 
4 



(26) 



1 rP-V 1 



standard terms 



new terms 
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where we distinguish between the terms coming from the standard GDW approach and the 
new contributions due to the presence of the zero mode. Finally the modified GDW heat 
trace expansion is obtained by plugging (fT9|) into (f20|) . 



M/3) 



^ Cn (^)/3" + Erf(^) 



Cn(K) 



J 

Jn 



dx Cfi ^ x : x ^ 



where the Seeley coefficients c n {K) are given by the spatial integral of the Seeley densities 
Cn(x,x). Therefore, the first three Seeley coefficients are 



co(K) 
c x {K) 

c 2 (K) 



I 



- (V(a?)) - Av , 

\ {V"W) + \ ((^(x)) 2 > + 4v (V(x)fi(x)) 



(27) 



standard terms 



new terms 



where again we point out the new terms introduced in this context. Here I is the length of 
the interval 0, which eventually we make to tend to infinity. In any case, subtraction of the 
-Ko-heat trace from the K-he&t trace expansion amounts to dropping the cq{K) coefficient, 
and we find: 



h K {P)-h Ka {P) 



^c n (iO/r-f +Erf( Uv / ^) . 



(28) 



n=l 



A physical application: computation of the one-loop kink 
mass shift 



3.1 Classical field theory models and kinks 

The notion of kink arises in the (l+l)-dimensional relativistic scalar field theory context 
[12\ [T31 [H] . In this framework the dynamics is governed by the action 



Sty] 



1 dip dtp 
2 dy^ dyv 



Uty(y»)] 



Here, ijj(y^) : R > — > R is a real scalar field; y° = r and y 1 = y are local coordinates 
in R 1 ' 1 , equipped with a metric tensor = diag(l, — 1), fjL,v = 0,1. We shall work in a 
system of units where the speed of light is set to one, c = 1, but we shall keep the Planck 
constant h explicit because we shall search for one-loop corrections, proportional to H, to 
the classical kink masses. In this system, the physical dimensions of fields and parameters 



arc: 



[h] = [S] = ML, [ Vtl ] = L, M 



U] = ML . In this type of models we can 



always identify two special parameters, rrid and 7,2, to be determined in each case, carrying 
the physical dimensions: [mj = L^ 1 and [7J = M~?L~2. 
non-dimensional coordinates, fields and potential: x^ 



u{4>) 



which allow us to introduce the 

m d y^, x = t, xi = x, 4> = 7dV, 



U(ip). The action and the "static" part of the energy are also proportional to 



dimensionless action and energy functionals, namely: 



1 



s\ 



1 

7 d 



dx°dx 1 



1 



2 dxu dx^ 



U[<t>{x»)\ 



(29) 
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%] = ^E[<f>] = ^ /dx 

"^d ^d 



lit) ^ 



(30) 



where we shall assume that U(4>) is a non-negative twice-differentiable function of 4>: U{4>) £ 
C 2 (R) and U{4>) > for G R. The configuration space C of the system is the set of field 
configurations C = {4>(to,x) G Maps(R 1 ,R)/E , [(/)] < +oo}. The static and homogeneous 
solutions of the field equations 



dt 2 



8P_ 

dx 2 



4>(t,x) 



5U 



(t,x) 



(31) 



correspond to the vacua of the model, M = {<f>® / U((f>^) = 0}. Topological solutions, 
kinks in this context, are defined as localized non-singular solutions of the field equation 
whose energy density, as well as being localized, has spacetime dependence of the form: 
e(t, x) = e(x — vt), where v is some velocity vector. The static kink solution 4>k{x) is a BPS 
solution connecting two vacua 00 and which satisfies the first-order ODE 



dx 



±^/2U{^ K ) . 



(32) 



Small perturbations around the kink and vacuum solutions are described by the spectral 
problem associated with the vacuum and kink fluctuation operators 



d 2 



+ 



d 2 u, 



(0i 



d 2 d2jJ ^ , m 

K = -— + ^[^{X)] 



dx 2 d(p 2 ' dx 2 

which adopt respectively the form ([5]) and (pQ) provided that we fix 



d 2 u, 



0i 



V(x) 



d 2 u 



4>k{x)] - v 2 



3.2 One-loop kink mass shift and the zeta function regularization 

Dashen, Hasslacher and Neveu (DHN) solved the problem of computing the shift in the kink 
mass induced by kink fluctuations in the one-loop order of the sine-Gordon and )i+i 
models in [T5] . The authors wrote the one- loop kink mass shift as the sum of two contributions 
AE((()k) = A£?i(0jj-) + AE2(4>k)'- the kink Casimir energy (mode-by-mode subtraction of 
the zero point vacuum energy) and the mass renormalization counter-term. In order to apply 
the DHN procedure effectively, it is necessary to know the eigenvalues of the bound states and 
the scattering wave phase shifts of the K operator. For a generic scalar field theory model 
we lack this spectral information. Alternatively a zeta function regularization can be used in 
the previous scheme, see [H |9l HH [T7] . The vacuum energy induced by quantum fluctuations 
is first regularized by assigning it the value of the spectral zeta function of the Xo-operator 
(a meromorphic function) at a regular point in s G C: 

A£(<^) = ^GcoH) -> AE&®)[s] = ^-f (4V Ck (s) , (33) 

where \x is a parameter of dimensions L~ 1 introduced to keep the dimensions of the regularized 
energy right. We stress that s = — \ is a pole of this function. The same rule is applied to 
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control the kink energy divergences by means of the spectral zeta function of K and therefore 
the kink Casimir energy AEi((J)k)[s] is regularized in the form: 



A = AE (<j> K )[s] - AE (^)[s] 



mi 



(Ck(s) - (ko(s)) , (34) 



On the other hand the energy due to the one- loop mass counter-term AE2(4>k)[s] can be 
also regularized by the zeta function procedure: 



AE 2 (<j> K )[s] = -^(V(x)) 



mi 



1- 

i-+oo i r(s) 



(k (s + 1) . 



(35) 



Finally, we write the zeta function regularized DHN formula: 

A%)= lim i AE 1 (ct> K )[s}+ lim x AE 2 ((/)k)[s] . (36) 

s— >— ^ s— >— ^ 

The Mellin transform expresses the relation (|34p in terms of the kink and vacuum fluctuation 
operator heat kernels: 



AE^k^s] 



2 



2\ s 



1 



dp P"- 1 (h K (P) - h Ko (J3)) 



(37) 



my T(s) 

which allows us to apply the modified GDW approach introduced in the previous section 
This leads to 



A£i(^)[s] 



h ll ( V 



+ 



mi 



E 

n=2 



1 {V{x))T[s + \] 



2 T[s + \] 



yfkx v 1+2s T[s] y/nv 2s T[s] 



+ 



c n (K) T[s + n 



2n+2s-l 



yjftv 2s sT[s] 



where we have used the explicit expression (|27p of c\ (K) to write the contribution of the first 
term of the series. Because the contribution of the regularized one-loop mass counter-term 
AE2(4>k)[s], see (|35|). can be written as 



m; 



(V(x))T[s + ±] 1 
V^F T[a] v 2 ^ 



^ oo 



(38) 



we end from (|36[) with the renormalized one-loop mass shift formula derived from the modified 
asymptotic series expansion of the if-heat function: 

AE^k) 

Either because of computational restrictions in the estimation of the Seeley coefficients or 
because of the asymptotic nature of the series, (|38p is truncated to a finite number Nt of 
terms, which provides us with 



7r 8ir 



n=2 



as a good estimation of the quantum kink mass correction. The most important advantage 
of this procedure is that there is no need for detailed information about the spectrum of K 
to calculate the mass shift. 
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3.3 An example: the sinh 4 model 



The potential term, which determines the dynamics in this model, is given by: 

U (</>) = - (sinh 2 </)- l) 2 . 

This function has two absolute minima located at cf^ 1 ' = — arcsinh 1 and 4>^ = arcsinhl, 
which play the role of vacua of the model. From the first order ODE (|32p we can identify the 
static kink solitary wave 

tanh x 

<Pk\x) = arctanh— -^=— (40) 

which connects the two vacuum points of the model. The fluctuations over this solution 
can be described by the spectral problems associated with the vacuum and kink Hessian 
operators, which are respectively: 

K = -^+A , K = -— r + 2 + 



dx 2 dx 2 (1 + sech 2 cc) 2 1 + sech 2 x 

The previous operators follow the form ([I]) and §5§ established in the first section, sim- 
ply taking v 2 = 4 and V(x) = — 2 ^^se^x^ ^ ' one ~l°°P m ass shift for the kink 
can be computed by the formula (j39|) where the zero mode fo(x), the normalized spatial 
derivative of the kink (|4"0j) . is fo(x) = , 4y/2 and the Seeley coefficients 

V3V2arccosh3-4(3+cosh(2a;) 

c n (K) are determined by means of the recurrence relations (|24p . The Mathematica pro- 
gram KinkMassQuantumCorrection_Modified.nb, which can be download at the web page 
http://campus.usal.es/~mpg/General/Mathematicatools, automatizes this calculation. The 
Seeley coefficients for the present case are listed in Table [TJ 



Table 1: Seeley coefficients and partial kink mass shift estimations 



n 


Cn(K) 


N t 


AE(<l> K ;N t )/hi* 


1 


7.47870 


1 




2 


7.82708 


2 


-0.714477 


3 


5.71901 


3 


-0.728699 


4 


3.15228 


4 


-0.732619 


5 


1.36104 


5 


-0.733888 


6 


0.482077 


6 


-0.734338 


7 


0.14436 


7 


-0.734506 


8 


0.0375685 


8 


-0.734572 



In this Table we have also specified the values of the one-loop mass shift obtained from the 
modified asymptotic series (|39p for different values of the truncation order, Nt. These data 
have been depicted in Figure 1. For sake of comparison the same data extracted from the 
standard asymptotic expansion are also included, see [9]. Notice that the modified method 
exhibits better convergence properties than the standard procedure. 
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Figure 1: One-loop kink mass shifts estimated by means of the modified and standard asymp- 
totic series for several truncation orders. 

In conclusion, the one-loop mass shift for the kink (|40p in the sinh 4 (f> model is given 
by the value —0.734572, which has been obtained from the modified heat trace expansion 
truncated at the order N% = 8. The asymptotic behaviors of the partial sums associated with 
the modified and standard GDW heat trace expansions are illustrated in Figure 2. Notice 
that the modified approach reproduces the correct analytical asymptotic behavior as opposed 
to the standard procedure. 

Figure 2: Graphics of the partial sums of the heat trace expansion from the modified (left) 
and standard (right) GDW approach. 
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